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1. Introduction 



String theory as a possible candidate for quantum gravity proposal implies existence 
of extra spacetime dimensions. In fact, existence of possible extra spatial dimensions 
has been a principal key to address the long-lived Hierarchy problem p]-[2]. The 
first attempt to adopt extra dimensions in theoretical physics was made by Kaluza 
and Klein to unify gravity with electromagnetism in a 5-dimensional spacetime model. 
Currently the idea of having extra dimensions has been widely accepted by theoretical 
physics community and there is also hope to detect some footprints of these extra 
dimensions in experiments such as the LHC After revolutionary report on current 
positively accelerating phase of the universe expansion [3], possible modification of 
General Relativity has been in the center of some important research programs. One 
may modify the geometric part of the Einstein's equations by allowing to have extra 
spatial dimensions in the framework of braneworld scenarios These attempts 

have been applied successfully to explain the issue of dark energy and interpretation 
of new astrophysical data. These are just parts of possible motivation to introduce 
extra dimensional scenarios in theoretical physics. Nevertheless, the presence of extra 
dimensions would remain concealed some how, may be because the extra dimensions 
are compactified to small radius with the size around Planck scale. Extra dimensional 
scenarios open also new perspectives, at least on phenomenological ground, to quantum 
gravity proposal. A common feature of all approaches to quantum gravity, such as 
string theory, noncommutative geometry, loop quantum gravity and Doubly Special 
Relativity, is modification of the standard dispersion relation (see for instance [5], [6], 
[7] and references therein). In fact, in the study of loop quantum gravity, Doubly 
Special Relativity and models based on non-commutative geometry there has been 
strong indication to modify the standard energy-momentum dispersion relation. In 
this framework issues such as possible violation of the Lorentz invariance in quantum 
gravity regime are studied too (see [8] and references therein). Recently modified 
dispersion relations are formulated in model universes with extra dimensions by focusing 
on black hole thermodynamics in these models [9] . Here we focus on another important 
issue: the block body radiation in a model universe with large extra dimensions and by 
incorporating possible quantum gravity effects encoded in modified dispersion relations. 
A black body is an ideal body which allows the whole of the incident radiation to pass 
into itself and absorbs within itself this whole incident radiation. One cannot ignore 
the significance of study of black body spectrum and its possible modification since 
it could be as a reliable source to give some information about the remnant of early 
universe radiation which could be detectable today as cosmic background microwaves 
(CMB). In fact, the cosmic microwave background radiation observed today is the 
most perfect black body radiation ever observed in nature, with a temperature of 
about 2.725 K. It is a snapshot of the radiation at the time of decoupling between 
matter and radiation in the early stages of the universe evolution. Therefore, any 
modification of black body spectrum due to quantum gravity effects can open new 
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windows to know more about quantum gravitational features of the early universe 
through study of CMB spectrum. This sort of investigation has been considered by some 
authors recently pU]. Also formulation of black body radiation in model universes with 
extra dimensions, but in the absence of quantum gravity effects, has been considered 
by some authors [11] . The crucial feature of our work is the combination of MDR 
(as a phenomenological outcome of quantum gravity proposal) with possible existence 
of extra dimensions in problem of black body radiation. This study is an attempt 
to take a small step toward a deeper understanding, at least phenomenologically, of 
quantum gravity in a model universe with large extra dimensions. In this framework we 
investigate the impact of modified dispersion relations (MDR) on black body spectrum 
in a model universe with large extra dimensions. We establish modified form of the 
Planck distribution, equipartition theorem, Stefan-Boltzmann law and Wien's law in 
this setup. Also modified Debye law and modified Dulong-Petit law are obtained and 
the corresponding results are interpreted. Then we study entropy and specific heat of 
black body radiation in this framework. The important ingredient and generic feature 
of these modified quantities is that these modifications are themselves temperature- 
dependent. So, it seems that quantum gravity effects encoded in modified dispersion 
relations are temperature-dependent for black body radiation and also related issues. 
This feature may be related to the foamy structure of space-time at the quantum gravity 
level, and taking the quantum structure of spacetime to be as a granular media which 
has some interaction (through quantum fluctuation of the background metric) with the 
events which are happened in it. So, these effects can be attributed to some yet not 
well-understood features of quantum spacetime at the Planck scale. We note also that 
the assumption of black body in a solid box, to some extent, helps us to know more 
about the modified terms within MDR in vibrating solids. Therefore the modified heat 
capacity of solids would be achieved through Dulong-petit model in high temperature. 
In which follows, we set h = ks = c = 1 for simplicity. 

2. Black body radiation in a model universe with extra dimensions 

In order to study a black body radiation and it's features in a D-dimensional space (D 
spatial dimensions), we consider a conducting cavity with the shape of a D-dimensional 
cube of side L. This cavity is filled with electromagnetic radiation which is in thermal 
equilibrium with walls at temperature T and is linked with the outside by a small hole. 
The electromagnetic radiation inside the cavity can be treated as standing waves. We 
are going to seek for the functional form of energy density for one-dimensional cavity 
that is placed in a D-dimensional universe (see Ref. pj] for a discussion about validity of 
this assumption). The most proper method for investigating the black body radiation 
is deriving the total number of modes in order to obtain the energy density. Let us 
write the components of the electric field in a cavity by choosing a system of orthogonal 
coordinates with origin at one of the cavity's vertices. To satisfy boundary conditions, 
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at the walls of the cavity the parallel components of the electric field must vanish. 
Taking these considerations into account, we have the following relation in terms of 
wave number 

TH = —L (1) 

7T 

where ki is the component of wave vector and rij = 1,2,... represents the possible modes 
of vibration. The total number of modes in one- dimensional space is given by 

N(ki)dki = -dh (2) 

7T 

which can be converted to the following form in a spacetime with D spatial dimensions 

V 
2^ L 

where V = L D is the volume of the cavity and dVk is the infinitesimal element of volume 
in the fc-space [TTJ. Coefficient (D — 1) is related to polarization states and the presence 
of 2 D in the denominator of equation (3) is due to the fact that we should consider only 
the positive part of the volume in fc-space. The D-dimensional element of volume in 
terms of v is given by 



N{k)dk = {D-l) 7 — 15 dV k (3) 



D 
71 2 



dV v = 2^{2?r) D v^dv (4) 

By using equation (4), now equation (3) in z/-space can be rewritten as follows 

N{v)du = 2{D-l)V^v^dv. (5) 

This relation expresses the total number of modes with frequencies between v and u+du. 

In the next section, by using the average energy per mode, we will construct the 
spectral energy density to explain black body radiation in the presence of quantum 
gravity effects in a spacetime with D spatial dimensions. 



3. MDR-modified black body radiation in a model universe with D spatial 
dimensions 

The modified dispersion relation in a model universe with D spatial dimensions can be 
written in the following form (see [5]- [9] for more details) 

f = f{E, m; L P ) ~ E 2 - fi 2 + aL 2 E 4 + a'L 4 p E 6 + 0(L 6 p E 8 ) , (6) 

where Lp, the Planck length in spacetime with D spatial dimensions, depends on 
spacetime dimensionality (see [9] for instance), and / is a function that gives the exact 
dispersion relation [5]. On the right hand side of equation(6) we have adopted a Taylor- 
series expansion for E The parameter a, is dimensionless and takes different 
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values depending on the details of quantum gravity candidates. The rest energy of the 
particle is m which is different from the mass parameter \i\ we expect fj, ^ m if the 
coefficients oti do not all vanish [5]. 

Based on the analysis of the previous section, the number of modes in the momentum 
interval between p and p + dp, is given by 



2vrf 



N(p)dp = (£>-!) V^p D ~ l dp. 



(7) 



To find the modified form of the total number of electromagnetic modes in MDR 
framework, we start with differentiating and applying a Taylor expansion to equation 
(6) to obtain 



dp = dE 



l + |«^ 2 +(-|« 2 + |«') L ^ 



(8) 



2 p V 8 2 

where we have saved only terms up to the forth order of the Planck length and the rest 
mass has been neglected. Substituting the expressions for p and dp, we find 



N(E)dE ={D- 1)V 



^E D - l (l + aLlE 2 + a'L A p E A 



X 



l+l«^ 2 +(-k+l«) L ^ 



dE 



(9) 



2 p V 8 2 

This is the modified number of electromagnetic modes in a cavity and in a model universe 
with D spatial dimensions. Note that the average energy per mode is independent on 
the dimensionality of space and it can be given by (see for instance [TTj ) 

v 



E 



er — 1 

Then, the modified spectral energy density is 



(10) 



p(v)dv — (D — 1 
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This is the modified Planck distribution. The influence of modified dispersion relation 
on spectral energy density in a model universe with large extra dimensions is obviously 
considerable. Figure 1 shows the variation of the spectral energy density versus v in 
Planck temperature for a model universe with D spatial dimensions but without MDR 
effect. This figure shows that when the number of spatial dimensions increases, the 
frequency at which the distribution maximizes (the location of the distribution peak) 
will attain a noticeable shift toward higher frequencies. Also the distribution becomes 
wider and its height decreases by increasing D. This shows that by increasing D, the 
number of modes located around the most probable frequency increases considerably. 
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Figure 2 depicts the variation of the modified spectral energy density versus v in D- 
spatial dimensions when quantum gravity effects are taken into account through MDR 
and temperature is fixed at the Planck temperature. While the location of peaks in 
this case are D-dependent in the same way as in Fig. 1, the height of picks increases 
by increasing D and the distribution width increases also with D. So, there is a sharp 
difference between the results obtained in the presence of quantum gravity effects and 
their absence. This feature may provide in principle a direct clue to test the underlying 
theory in the lab. 



p(v,T) 




v 



Figure 1. Spectral energy density as a function of frequency in a model universe with D spatial 
dimensions and without MDR effect. 

Comparison between these two figures shows that when quantum gravitational 
effects are taken into account, the frequency at which the distribution maximizes will 
attain a noticeable shift toward high frequencies. This is because of the impact of 
quantum gravitational terms on spectral energy density. Now we explain an important 
achievement of our analysis through figures 3 and 4. These two figures show the spectral 
energy density as a function of frequency with and without incorporation of quantum 
gravity effects for D = 3 and 8 (just for example) respectively. By incorporation of 
quantum gravity effects via MDR, in each figure the frequency at which the distribution 
maximizes has significant shift toward the higher frequencies and the height of the 
distribution increases too. Comparison among these two figures shows that as the 
number of spatial dimensions increases, the difference between two peaks (frequency 
shift) increases. Based on these figures, the effects of quantum gravity corrections get 
enhances as the number of spatial dimension's of the universe increases. 
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Figure 2. Spectral energy density as a function of frequency in a model universe with D spatial 
dimensions and in the presence of quantum gravity effect encoded in MDR. 
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Figure 3. Spectral energy density as a function of frequency for D=3 in two different regimes. 

In 4- dimensional space-time, classical equipartition theorem states that average 
energy per electromagnetic modes of black-body radiation is T (with ks = 1). In 
ordinary quantum mechanics, this average is given by f . Now in a model universe 
with large extra dimensions and in the presence of quantum gravitational effects encoded 
in MDR, this average generalizes to 
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Figure 4. Spectral energy density as a function of frequency for D=8 in two different regimes. 



2 -5 3(D-1), „5 (D-l ; >. , , , 
a 2 (— + -- — - — -) + a'{- + - — - — -) )L> 4 



(12) 



3.1. Modified Rayleigh- Jeans law 



The classical Rayleigh- Jeans law describes the spectral energy density of electromagnetic 
waves in black body radiation, which fits with experimental results at low frequencies, 
while in high frequencies it tends to infinity known historically as the Ultraviolet 
Catastrophe. The modified form of the Rayleigh- Jeans law in our framework is written 

as 



5 3 (Q-1K ,/5 (g-i; 

a + a — 

v 8 2 2 /^V2^ 2 



L> 4 T. 



(13) 



Thus, the modified Jean's number can be written as 



r(f) 



3 (D-l)x r2 2 
2 + 2^) aL ^ + 



Q 



,/-5 



2 ,-5 , 3CD-l)x ( 6 

2 2 / 2 



4,. 4 



(14) 



Assuming a, are positive quantities, the Jeans's number as number density of modes in 
a given frequency, increases in the presence of quantum gravity effects. We note that 
since the Reyleigh- Jeans theory has a classical basis, we have used ksT (with ks = 1) 
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as mean energy per mode in the above calculations. Nevertheless, in a general frame- 
work we should take into account the modification of equipartition theorem and use 
equation(ll) as our starting point. The result is the same as equation (14). 



3.2. Modified Stefan- Boltzmann law 

In ordinary quantum mechanics with D = 3, the total energy density of black body 
radiation is proportional to the forth power of temperature, that is, p(T) = aT A . 
In a higher dimensional world, when ordinary quantum mechanics is taken into account, 
the total energy density is p(T) = a£>T D+1 [IT] . Now, in a model universe with large 
extra dimensions and in the presence of quantum gravity effects encoded in the MDR, 
we integrate equation (11) over frequency to find 



pen 



[D - 1)2At d+1 



r 



3 CD-I) 
2 2 



+ 
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C(D + i)r(D + i) 

5 3 (D — 1)\ i ,(h i D — I 



8 2 



+ «'(- + 



aL p 2 T 2 ((D + 3)T(D + 3) 
L p 4 T 4 C(D + 5)r(D + 5)l(15) 



In comparison with p(T) = odT d+1 , the coefficient old in the presence of quantum 
gravity effects is given by 



a D 



(D - 1)2tt 



C(D + 1)T(D + 1) 



(D-l) 



2 + 2 



+ 



5 3 (D-i; 
8 2 2 



,/5 D-l 

+ a (-H 

V2 ^ 2 



aL p 2 T 2 ((D + 3)r(D + 3) 



L p 4 T 4 C(D + 5)r(D + 5) (16) 



where C and T are the Riemann Zeta function and gamma function respectively. The 
Radiancy is defined by the total energy emitted by a black body per unit area per 
unit time. The energy density p(T) is proportional to the Radiancy R(T) by a pure 
geometric factor [11] 



D ~ 



R(T) 



1 



r^±i)2 v ^F 



(17) 



Therefore, the modified Stefan-Boltzmann law may be expressed as 

|aL p 2 T 2 C( J D + 3)r( J D + 3) 



R( T ) = ^^ m ^-T D+1 \aD+l)T(D + l) + 



+ 



3 (D-l) 
2 2 



a- - 
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L p 4 T 4 C(D + 5)r(D + 5) (18) 
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Now, o~m as the modified Stefan-Boltzmann factor in a model universe with D spatial 
dimensions and in the presence of quantum gravity effects encoded in MDR can be 
written as 



7T 2 



aL p 2 T 2 ((D + 3)T(D + 3) 



+ 



a 



5 3 ( /; i ) 

2 2 > '2 2 



)+o/(l + °-±) L/T 4 C(^ + 5)r(D + 5)kl9) 



We note that the modified Stefan-Boltzmann factor in this framework is not constant 
and varies with temperature. Also the correction terms contain only the even power 
of temperature. This temperature dependence of modifications is a generic feature and 
probably has its origin in the quantum fluctuation of background metric in quantum 
gravity regime. 



3.3. Modified Wien's law 

In ordinary quantum mechanics and in three space dimensions, there is an inverse 
relationship between the wavelength of the peak of the emission distribution of a 
black body and it's temperature, which is called Wien's displacement law. This law is 
expressed as \ max = |;, where b is Wien's constant and T is the absolute temperature of 
the black body when X ma x is the distribution peak wavelength. In a higher dimensional 
universe and by incorporating the quantum gravity effects, we should rewrite the 
modified energy density, equation (11), as a function of wavelength. Then we find 

2ttt 1 [ ,3 (D-1)n ,,1 r o/-5 3 



/ m x / ^ ,Z7T2 1 /O I U i J \ 9 1 9 / — O 



+ 0i + — 5 — ) 




(20) 



In order to find the wavelength of maximum emission, we differentiate p\(T) (expressed 
up to the second order of Lp) with respect to A, and then by using the approximation 
e^r ps 1 + and putting the derivative equal to zero, we find 

1 (36BA - 108FA 2 - 8 + 12y/3y/AB 3 A - B 2 - 18BAF + TIF 2 A 2 + 4FA) * 

- 2 - WA ~ l r (21) 

3 A(36BA - 108FA 2 - 8 + UVsVW^A - B 2 - 18BAF + TIF 2 A 2 + AFA) 3 

where A, I? and F are defined as 
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A=-T(D+1), B = -T(D + 3)(^+ i -^ Y ^-)aL 2 P , 



Note that the first term in (21) gives the standard Wien's law in the absence of 
quantum gravity effects. The extra terms in (21) are quantum gravity corrections and 
as usual these corrections are temperature dependent. 



4. Entropy and specific heat 

Now we calculate entropy and specific heat of black body radiation in a model universe 
with D spatial dimensions and in the presence of quantum gravity effects encoded in 
MDR. We start with total energy in cavity as 



U(T) = Vp(T) = V— ^ { C(D + l)r(£> + 1) 



Q 



+ (| + )aL/T\(D + 3)T(D + 3)+ 

5 3(D-1 
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Then the specific heat and entropy of black body radiation are expressed respectively 
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.(24) 
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Figures 5 and 6 demonstrate the variation of specific heat versus temperature 
without and with MDR respectively. As figure 5 shows, the heat capacity decreases 
by increasing the number of spatial dimensions. Comparison between figures 5 and 6 
shows that for a given temperature, the value of specific heat with quantum gravity 
effects is larger than the corresponding value in the absence of these effects. In other 
words, incorporation of quantum gravity effects through MDR increases specific heat of 
black body radiation for a given temperature. 




T 

Figure 5. Black body radiation specific heat in a model universe with D spatial dimensions without 
MDR effect. 




Figure 6. Black body radiation specific heat in a model universe with D spatial dimensions in the 
presence of MDR effect. 
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Figures 7 and 8 demonstrate the behavior of black body entropy versus temperature 
for different numbers of spatial dimensions without and with MDR effects respectively. 
Figure 7 shows that the entropy decreases as the number of spatial dimensions increases. 
Comparison between figures 7 and 8 shows that quantum gravity effects increase the 
entropy content of black body radiation in a given temperature. 




Figure 7. Black body entropy in a model universe with D spatial dimensions without MDR effect. 




T 

Figure 8. Black body entropy in a model universe with D spatial dimensions in the presence of 
MDR effect. 
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5. Modified Debye's law 



Imagine a black body as a solid box, which is compound of vibrating atoms. These 
vibrating atoms are recognized as the so-called phonons in this box. The method 
adopted in this section resembles more or less calculation of number of electromagnetic 
modes in black body radiation. Nevertheless, there are some discrepancies between these 
two cases since in a vibrating solid there are three types of waves (one longitudinal and 
two transverse) in three spatial dimensions while there are only two polarization states 
for electromagnetic modes. So, if we generalize this situation to a model universe with 
D spatial dimensions, we would have D possible polarization states for phonons. Unlike 
photon, phonon frequency is bounded by upper limit. In fact, the total vibrational 
states for phonons is D x iV where N is the number of atoms in vibrating solid. Since 
there is no limitation on the number of phonons in each energy levels, phonon could 
be considered as bosons. Therefore, in a model universe with D spatial dimensions, we 
have 

U{y) = DL D -^ I -rdu (25) 

r( Y )Jo (eT-l) 

which determines the total energy of vibrating solid in a model universe with D spatial 
dimensions. By setting x — ^, we obtain 

/ T \D « m „ x d 
U{X) = D >TN(-) I — & (26) 

where Tp is introduced as Debye temperature which is equal to 

Now we incorporate quantum gravity effects through MDR. Based on our previous 
analysis, we have 



/ T \ D r x 




This relation leads us to the modified Debye law straightforwardly. Since the integral 
cannot be calculated analytically, we consider some limiting cases. The temperature of 
a Debye solid is said to be low if T T D . In this limit we find 



U(T) = D 2 TN^—^j |C(£> + l)r(D + 1) + (- + L_J.y Lp i T 2 ((D + 3)F(D + 3) 
2 ( - | + +«'(! + W + 5)rp + 5) } .(28) 
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By differentiating this relation with respect to T, we find 
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This is the modified Debye law at very low temperature in the MDR framework. Note 
that the first term in the right hand side gives the standard Debye law at very low 
temperature in the absence of quantum gravity effects. It's worth mentioning to note 
that modified specific heat changes with temperature. So, quantum gravity corrections 
in this case are temperature dependent as usual. Figure 9 shows the heat capacity of 
solid based on Debye model at low temperature as a function of temperature. Based on 
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Figure 9. Debye specific heat of solid at very low temperature in two different regime for D=3 and 
at a fixed To . 



Debye model, the value of specific heat for solid at low temperature with MDR is larger 
than the corresponding value in the absence of MDR at a given temperature. 
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5.1. Modified Dulong-petit law 

Dulong and Petit found that the heat capacity of a mole of many solid substances 
is about 3R, where R is the universal gas constant. An equivalent statement of the 
Dulong- Petit law in modern terms is that, regardless of the nature of the substance 
or crystal, the specific heat capacity C of a solid substance (measured in Joule per 
Kelvin per Kilogram) is equal to ^ where M is the molar mass (measured in Kilogram 
per Mole). Thus, the heat capacity per mole of many solids is 3R. Actually this 
statement is valid only at high temperature and in this case the specific heat is assumed 
to be temperature-independent in classical view point. Now we reconsider this model 
in the MDR framework in order to see how quantum gravity effects modify the classical 
Dulong-Petit law. At high temperature when T > Tq, we can use the approximation 
e x — 1 PS x. Therefore, based on Eq. (27), the total energy of vibrating solid is given by 



Um . DTN + ^4 + ^±) aLr + + 1^. 

By differentiating with respect to T, we find 

Cv = DN+ \JdTT) NTd {2 + ^- ) aL p + 2 [ a 1^ + 2^—) 

This is the modified Dulong-petit law. In the absence of quantum gravity effects, 
ordinary quantum mechanics result will be achieved, which states that the heat capacity 
is independent of temperature. However, in quantum gravity domain one cannot ignore 
temperature dependence of heat capacity as is evident in correction terms. These 
temperature-dependent corrections reveal some phenomenological aspects of quantum 
gravity proposal. These features may provide also a direct basis to test quantum gravity 
in the lab. Figure 10 shows the heat capacity of solid based on Dulong-petit model at 
high temperature clS du function of temperature. 

Specific heat of solid in Dulong-petit law within MDR has linear dependence on 
temperature (at least up to our adopted approximations), while without MDR it is 
temperature-independent. 



Black body radiation in a model universe with large extra dimensions and quantum gravity effectsYI 



o- 

0.01 0.02 0.03 0.04 0.05 0.06 

T 

Figure 10. Dulong-petit specific heat of solid in two different regimes for D=3 at a fixed Tr>. 

6. Conclusion 

It is now a well-established fact that the Cosmic Background Microwave Radiation has a 
perfect black body spectrum. Modification of this spectrum due to tiny quantum gravity 
effects and possible detection of these footprints in CMB spectrum via observation opens 
new windows to know more about quantum gravitational features of the early universe 
and some phenomenological aspects of quantum gravity itself. Also it is interesting in 
essence to know the black body spectrum if our observable universe could be braneworld 
in nature. Based on these motivations in this paper we studied the problem of black 
body radiation in a model universe with extra dimensions by incorporating quantum 
gravity effects encoded in modified dispersion relations. We calculated with details 
some important physical and thermodynamical quantities attributed to black body 
radiation. As a generic result we have shown that quantum gravity corrections are 
generally temperature dependent. We have shown that in the absence of quantum 
gravity effects, when the number of spatial dimensions increases, the frequency at which 
the modified Planck distribution maximizes (the location of the distribution peak) will 
attain a noticeable shift toward higher frequencies. This distribution becomes wider 
and its height decreases by increasing the number of spatial dimensions. This means 
that by increasing D, the number of modes located around the most probable frequency 
increases considerably. When quantum gravity effects are taken into account through 
MDR the location of peaks are D-dependent in the same way as in the absence of 
quantum gravity effects, but the height of picks and the distribution width get increase 
by increasing the number of spatial dimensions. This feature shows that there is crucial 
difference between the results obtained in the presence of quantum gravity effects and 




Without MDR 
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their absence and this may provide, in principle, a direct clue to test the underlying 
quantum gravity candidate. We have shown also that the effects of quantum gravity 
corrections get enhances as the number of spatial dimension's of the universe increases. 
As another important result, we have shown that the Jeans's number as number density 
of modes in a given frequency increases in the presence of quantum gravity effects. The 
modified Stefan-Boltzmann factor in this framework is no longer a constant and varies 
with temperature. The correction terms to Stefan-Boltzmann law in this framework 
contain only the even power of temperature. We obtained the modified Wien's law which 
shows a small, temperature-dependent shift in X max . In the next step, we have shown 
that incorporation of quantum gravity effects through MDR increases the specific heat 
of black body radiation for a given temperature. The entropy of the system decreases as 
the number of spatial dimensions increases. However, quantum gravity effects generally 
increase the entropy content of black body radiation in a given temperature. Based 
on the modified Debye model presented here, the value of specific heat for solid at low 
temperature with MDR is larger than the corresponding value in the absence of MDR 
at a given temperature. We obtained also the modified Dulong-Petit law. While in 
ordinary quantum mechanics picture of this law the heat capacity is independent of 
temperature, in quantum gravity domain one cannot ignore temperature dependence of 
heat capacity as is evident in correction terms. Once again this temperature dependent 
feature can be tested essentially in future experiments. Finally we have shown that up 
to adopted approximations, the specific heat of solid in Dulong-petit law within MDR 
framework has linear dependence on temperature but without MDR it is temperature- 
independent. 
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